Orbital order in Lao.sSri 5 Mn0 4 : beyond a common local Jahn- Teller picture 



Hua Wu, 12 f]C. F. Chang, 1 O. Schumann, 1 Z. Hu, 1 ' 3 J. C. Cezar, 4 T. Burnus, 1 
N. Hollmann, 1 N. B. Brookes, 4 A. Tanaka, 5 M. Braden, 1 L. H. Tjeng, 1 ' 3 and D. I. Khomskii 1 

II. Physikalisches Institut, Universitdt zu Koln, Zillpicher Str. 77, 50937 Koln, Germany 
2 Department of Physics, Fudan University, Shanghai 200433, China 
3 Max Planck Institute for Chemical Physics of Solids, Nothnitzerstr. 40, 01187 Dresden, Germany 

^European Synchrotron Radiation Facility, Boite Postale 220, 38043 Grenoble Cedex, France 
5 Department of Quantum Matter, ADSM, Hiroshima University, Higashi- Hiroshima 739-8530, Japan 

The standard way to find the orbital occupation of Jahn- Teller (JT) ions is to use structural 
data, with the assumption of a one-to-one correspondence between the orbital occupation and the 
associated JT distortion, e.g. in 06 octahedron. We show, however, that this approach in principle 
does not work for layered systems. Specifically, using the layered manganite Lao.sSri.sMnC^ as an 
example, we found from our x-ray absorption measurements and theoretical calculations, that the 
type of orbital ordering strongly contradicts the standard local distortion approach for the Mn 3+ 06 
octahedra, and that the generally ignored long-range crystal field effect and anisotropic hopping 
integrals are actually crucial to determine the orbital occupation. Our findings may open a pathway 
to control of the orbital state in multilayer systems and thus of their physical properties. 
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Rich physical properties of 3d transition metal (TM) 
oxides are largely related to an interplay among charge, 
orbital, spin and lattice degrees of freedom [TJ[2]. In par- 
ticular, the orbital pattern plays a crucial role in mag- 
netic exchange [3 and in many other properties. There- 
fore large attention is nowadays devoted to establishing 
the detailed type of orbital order (00) in oxides having 
degenerate orbitals, and different methods can be ap- 
plied. The oldest but simple and still the most widely 
used one is based on the assumption that there is a 
one-to-one correspondence between the orbital occupa- 
tion and the local distortion of the nearest-neighboring 
anion cage around the TM ion, i.e. the local Jahn- Teller 
(JT) distortion. 

For the doubly-degenerate case of e g electrons one can 
characterize orbital occupation by the "mixing angle" 
|<9) = cos(<9/2)|3z 2 - r 2 ) + sm(6/2)\x 2 - y 2 ). Local de- 
formations of the oxygen octahedron around the TM ion 
can also be presented as a superposition of the tetrag- 
onal distortion Q% and the orthorhombic distortion Q2, 
so that the general deformation of the octahedra is 
|<9') = cos((9 / )|Q 3 ) + s'm(0 f )\Q 2 ). When using the JT 
mechanism of electron-lattice interaction, the standard 
assumption, always made, is that the "mixing angles" 6 
and 6' are the same, and that the orbital mixing angle 6 

is determined by the relation [3] tan(0) = ^3^, where 
/, m and s are the long, middle and short distances, re- 
spectively, from TM to the nearest-neighboring ligands, 
all representing a local JT distortion. This rule is usually 
true for isolated JT centers, and for cubic oxide materials 
where an additional superexchange mechanism may also 
get involved in determining an orbital state [TJ [4] . Even 
in the latter case, however, after such an orbital state 
is established, a local distortion will in principle adjust 
to the orbital occupation, such that in effect again the 



local distortion and the orbital occupation would be the 
same. But, as we will show below, this rule breaks down 
for certain layered systems, because not only local distor- 
tions but also the anisotropic long-range contribution to 
the crystal field and the electronic kinetic (band) energy 
play crucial role in determining the orbital state. 

The half-doped single-layer manganite Lao.5Sri. 5 Mn04 
turns out to be a crucial testing and battling ground for 
the modeling of the relationship between crystal struc- 
ture, local electronic structure and the magnetic prop- 
erties. On the basis of crystal structure data, resonant 
scattering measurements and theoretical calculations [5]- 
HO] , the claim was made that the orbital pattern in the 
charge-ordered (CO) state of this material with the CE- 
type antiferromagnetic structure involves the cross-like 
x 2 -z 2 /y 2 -z 2 type orbitals for the Mn 3+ ions. The re- 
sults seem to fit to the standard local distortion model: 
local structural data show that the 06 octahedra around 
the Mn 3+ sites are not elongated, but rather locally com- 
pressed, with four long and two short Mn-0 distances, 
thus 'justifying' the cross-like type of orbital occupation. 
Yet this finding is quite surprising in view of the fact 
that earlier various theoretical studies [TTHT5] , not know- 
ing the local structure around the Mn ions, proposed the 
rod-like 3x 2 -r 2 /3y 2 -r 2 type to explain the experimen- 
tally observed magnetic structure. Also an alternative 
interpretation [16] of the x-ray scattering data was not 
considered to be persuasive. 

We reinvestigated this question on high-quality single 
crystals of Lao.5Sri.5Mn04 [T7H2Q] , using polarization de- 
pendent soft x-ray absorption spectroscopy (XAS) at the 
Mn I/ 2 ,3 edges [21] and multiplet cluster calculations [22]- 
[24] . We also performed density functional calculations 
within the local-spin-density approximation (LSDA) and 
LSDA plus Hubbard U (LSDA+/7) [25] [26] with a full 
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FIG. 1: (Color online) Polarization-dependent Mn-L2,3 XAS 
spectra of Lao.sSri.sMnCU for E||c (red curve) and EJ_c 
(black curve) taken at 150 K. (b) The corresponding linear 
dichroism (LD) spectrum. Theoretical LD calculated for or- 
bital scenarios with (c) Mn 3+ 3x 2 -r 2 /3y 2 -r 2 or (d) Mn 3+ 
x 2 -z 2 /y 2 -z 2 orbital occupation, each with the Mn 4+ t 3 ^. 

structural optimization [27] , thereby staying close to the 
recently available structural data [9j[10]. All the exper- 
imental and computational details are given in the Sup- 
plemental Material. 

Figure[I](a) shows our polarization dependent Mn-L2,3 
XAS spectra of Lao.sSri.sMnC^ for E||c (red curve) and 
E±c (black curve), where E is the polarization vector of 
the light and c the crystallographic axis perpendicular 
to the ab-plane. The spectra were taken at 150 K, i.e. 
below the orbital ordering temperature Too of 217 K. 
Curve (b) depicts the corresponding linear dichroic (LD) 
spectrum, defined as the difference between the two po- 
larizations, i.e. XAS(E±c) - XAS(E||c). We note that 
our spectra are very different from those reported ear- 
lier [5j |28] , and we will show below that our spectra are 
truly representative for the Lao.5Sri. 5 Mn04 system. 

To obtain some basic information concerning the or- 
bital occupation of the Mn 3+ ions, we can make use of 
a sum rule on the LD spectrum [29 . Integrating the 
LD spectrum throughout the £2,3 range, we find read- 
ily a negative net value. Since the Mn 4+ ions have the 
half-filled t\ g shell configuration, their contribution to 
the LD integral is zero. The non-zero value of the in- 
tegral is then due to the Mn 3+ ions only. The negative 
value directly indicates that the e g holes have a more 
out-of-plane than in-plane character, meaning that the 
e g electrons have their charge density more in-plane than 
out-of-plane. This then favors the rod-like 3x 2 -r 2 /3y 2 - 
r 2 type of orbital occupation and directly rules out the 
x 2 -z 2 /y 2 -z 2 scenario. 

In order to further confirm the above experimental 
finding, we have simulated the Mn-L2,3 LD spectra using 



the well-proven configuration interaction cluster model 
which includes the full atomic multiplet theory [22H24] . 
Figure [TJc) shows the calculated LD spectrum for the 
Mn 3+ ions having the rod-like 3x 2 -r 2 /3y 2 -r 2 orbital oc- 
cupation while the Mn 4+ ions are kept in t\ g configura- 
tion. We can see that all features of the experimental 
LD spectrum are nicely reproduced, including the neg- 
ative value for the integral. As a check we also include 
curve (d) which depicts the LD spectrum for the Mn 3+ 
ions with the cross-like x 2 -z 2 jy 2 -z 2 orbital occupation. 
We notice significant discrepancies with the experimental 
spectrum, especially at the L 2 edge where the sign of the 
simulated LD is opposite to that of the experimental one. 
Obviously, the integral value of this simulation (positive) 
has also the wrong sign. 

We would like to remark that the excellent agreement 
between the simulation in curve (c) and the experimen- 
tal LD spectrum can be taken as evidence that we have 
been able to obtain spectra which are representative for 
the Lao.5Sri.5Mn04 system. We have taken care that the 
high-quality single crystal was cleaved in-situ to obtain 
well ordered and clean sample area, that charging prob- 
lems were avoided by using a small sample area, and that 
the polarization of the light rather than the sample was 
rotated in order to keep the same sample area to be mea- 
sured for optimal comparison between the spectra taken 
with the two polarizations. 

Thus an important aspect that emerges directly from 
our experimental LD and multiplet cluster calculations 
is that the Mn 3+ ions have the rod-like 3x 2 -r 2 /3y 2 -r 2 
orbital symmetry. This is contrary to a prediction based 
on the local JT distortion: from the very similar in-plane 
2.01 A (x) and the out-of-plane 1.99 A (z) Mn-0 dis- 
tances, and the shorter in-plane 1.90 A (y) [27 , one 
could have expected the x 2 -z 2 level of e g crystal-field 
to be lower. This implies a mechanism for stabilizing the 
rod-like 3x 2 -r 2 /3y 2 -r 2 OO in this layered manganite dif- 
ferent from the conventional local JT mechanism. 

To understand better this mechanism, we have studied 
the electronic and magnetic structure of Lao.5Sri.5Mn04 
and particularly the origin of 3x 2 -r 2 /3y 2 -r 2 OO, us- 
ing LSDA and LSDA+Z7 calculations [251 [26] within the 
space group Bmbm which allows for the occurrence of the 
experimentally observed CO and OO. Figure |2|a) shows 
the orbitally resolved partial density of states (PDOS) for 
the two inequivalent Mn sites, calculated using LSDA+/7 
with U =5 eV and a Hund exchange energy of 0.9 eV. An 
insulating gap of about 1.0 eV lies between the split e g 
levels. For the Mn 3+ site, a pure 3x 2 -r 2 level appears just 
below the Fermi level. The electron count is less than 1, 
i.e. 0.41e, due to the bonding with the oxygens, as re- 
vealed by the mixed states below -4 eV. For the Mn 4+ 
site, the amount of occupied e g states is only a little less 
than for the Mn 3+ , but it is highly mixed x 2 -y 2 and 
3z 2 -r 2 hardly with any orbital polarization. Note that 
in spite of the e^-electron derealization, Lao.5Sri.5Mn04 
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FIG. 2: (Color online) (a) PDOS of the Mn 3+ and Mn 4+ ions in the CE-type antiferromagnetic state calculated by LSDA+(7. 
(b) Charge density contour plot (0.05-0.4 e/A 3 ) of the Mn 3+ -0-Mn 4+ network in the xy plane for the occupied e 9 states within 
1 eV below Fermi level. The ferromagnetic zigzag chain is marked by the dashed line, (c) Charge density contour plot of the 
Mn 3+ -0-Mn 4+ chain in the xz and yz planes, (b) and (c) clearly show the occupied 3x 2 -r 2 /3y 2 -r 2 orbitals of the Mn 3+ ions. 



can still be well categorized into a site-centered CO-00 
system; our calculations show that this state is more sta- 
ble than a bond-centered Mn-O-Mn polaronic state, by 
70 meV/Mn. This site-centered CO-00 state for the 
half-doped case has also been confirmed by previous ex- 
perimental [19] and theoretical [30] studies. 

The difference between the Mn 3+ and Mn 4+ sites can 
also be seen in Figs. [2|b) and [2|c), where we show 
the e g charge density contour plot on the xy plane as 
well as on the xz and yz planes of the Mn-0 network. 
One can clearly observe the 3x 2 -r 2 /3y 2 -r 2 OO at the 
Mn 3+ sites. The Mn 4+ ions, on the other hand, have a 
nearly isotropic e g charge density distribution. The CE- 
type ground state magnetic structure can be understood 
via the Goodenough-Kanamori- Anderson superexchange 
mechanisms, and this is also confirmed by our calcula- 
tions: the calculations of total-energy difference between 
various magnetic structures allow us to derive that the 
ferromagnetic coupling within each zigzag chain is about 
-30 meV, and that the inter- zigzag- chain antiferromag- 
netic coupling is much weaker, being only about 2 meV. 
This agrees qualitatively with the inelastic neutron scat- 
tering study [15] . 

It is important to note that the above results remain 
qualitatively unchanged within a wide range of U values: 
0-8 eV. The LSDA (U=0) calculation gives an insulating 
gap of 0.5 eV for the CE-type magnetic structure and 
0.32e (0.05e) polarization for the Mn 3+ (Mn 4+ ) e g or- 
bitals. Switching on the J7, the gap increases up to 1.3 
eV for U=S eV and the corresponding e g orbital polar- 
ization is 0.50e (0.05e) for the Mn 3+ (Mn 4+ ). The 3x 2 - 
r 2 /3y 2 -r 2 type of OO is found to be [/-independent. 

As mentioned above, the obtained orbital occupation 
contradicts a common local JT picture. Two factors con- 
tribute to this. First, there are longer-range contribu- 
tions to crystal field, cf. e.g. Refs [31j[32]. Especially 
in this layered manganite the coordinations of the Mn 3+ 
ion for further in-plane neighbors are very different from 
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FIG. 3: (Color online) (a) A sketch of the surrounding of a 
Mn 3+ ion in Lao.sSri.5Mn04 illustrates a different contribu- 
tion of the further neighbors to the crystal field. The local 
Mn 3+ -0 bondlengths are also shown, (b) The Mn 3+ e g PDOS 
calculated by LSDA, projected onto the (3x 2 -r 2 ,y 2 -z 2 ) basis 
giving eigenstates (upper panel), and onto the (x 2 -z 2 ,3y 2 -r 2 ) 
basis, with mixed occupation of both orbitals (lower panel). 



the out-of-plane ones. The further in-plane neighbors 
are the Mn 4+ ions, which have higher valencies and are 
at shorter distances than the further out-of-plane neigh- 
bors, being the (Sr 2+ , La 3+ ) ions, see Fig. [3|a). This 
has important consequences for the competition between 
the rod-like 3x 2 -r 2 and the cross-like x 2 -z 2 to become 
the lowest in energy. Our LSDA calculations reveal that 
the rod-like 3x 2 -r 2 orbital lies 90 meV lower than the 
cross-like x 2 -z 2 orbital. As a consequence, the former 
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FIG. 4: (Color online) Variation of the Mn 3+ e\ orbital state 
with respect to the c-axis Mn-0 bondlength R z , in terms of 
the orbital mixing angle calculated by LSDA+[7 (U=5 eV). 
The error bar is ±4° for U = 0-8 eV. The insets show the defi- 
nition of the orbital states, and the real Mn 3+ -0 bondlengths 
with R z — 1.99 A at which the 00 state is practically the 
pure 3x 2 -r 2 /3y 2 -r 2 type as discussed in the main text. 



becomes an occupied eigen-orbital [33 , while the latter 
is not, i.e. the occupied e g state of the Mn 3+ ion is a mix- 
ture of x 2 -z 2 and 3y 2 -r 2 if one chooses these orbitals as 
basis, see Fig. [3jb). A possible second factor is that the 
rod-like orbitals are more in-plane, thereby maximizing 
the gain in kinetic energy in this layered system. This 
factor is especially important for the systems close to 
the localized-itinerant crossover, to which our system be- 
longs. The importance of these two factors is confirmed 
by recent model calculations [34] . 

Nevertheless, if the c-axis Mn 3+ -0 bondlength were 
increased still further, the Mn 3+ e x g state would start to 
deviate from the 3x 2 -r 2 state, as seen in Fig.|H If that 
bondlength exceeded 2.2 A, the state would become 
more like x 2 -z 2 . Note, however, that the actual stretch- 
ing of the Mn06 octahedra in Lao.5Sri.5Mn04 along the 
c-direction, 1.99 A , caused by the layered structure itself, 
is much less than this 2.2 A value [27] . 

To summarize, using the linear dichroism and cluster 
model simulations of soft x-ray absorption spectroscopy, 
and density functional calculations, we demonstrate that 
Lao.5Sri.5Mn04 displays the rod-like 3x 2 -r 2 /3y 2 -r 2 or- 
bital ordering at Mn 3+ sites. Contrary to a common local 
JT picture, the observed orbital occupation strongly de- 
viates from the one expected from the local distortion 
of Mn06 octahedra, i.e. the assumption that the orbital 
and distortion mixing angles are the same, = 0' , is vio- 
lated. We explain this by considering the contribution of 
further neighbors to the crystal-field splitting, and by the 
important role of the electron kinetic energy in a system 
close to an insulator-metal transition. We expect that 
our findings are not only applicable for the '214' oxides, 
but also for layered materials in general, including the 
highly topical artificial multilayer oxide systems. 

This work is funded by the DFG via SFB 608. 
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